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ABSTRACT 

Using arguments from two dimensional Yang-Mills theory and the collective 
coordinate formulation of the Calogero-Sutherland model, we conjecture the dy- 
namical density correlation function for coupling I and 1/7, where I is an integer. 
We present overwhelming evidence that the conjecture is indeed correct. 
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There has been much recent interest in the Calogero-Sutherland Model (CSM) 
of interacting fermions in one dimension [1,2]. One reason for the interest is that 
these models are related to quantum spin chains, with long range interactions 
between the spins [3]. The CSM models also have a close relation to random matrix 
theory[2,4] and fractional statistics[5,6,7]. 

The Hamiltonian of the CSM on the circle is given by 

The spectrum of H has been long known [2], but until recently the dynamical cor- 
relation functions remained unknown. The basic problem is that while in principle 
it is possible to calculate the eigenf unctions for the CSM [2], in practice it is a com- 
putational nightmare. A big break in this direction came last year, when Simons, 
Lee and Altshuler calculated the dynamical two point density correlation function 
for three special values of the coupling. They exploited the fact that for these val- 
ues of the coupling, the CSM was related to random matrix theory of orthogonal, 
unitary or symplectic matrices. 

In this letter we will generalize the results of Simons Lee and Altshuler for CSM 
couplings I and 1//, where / is any integer. Our analysis is based on some straight- 
forward observations about the form of the measure in the correlation functions. 
Using arguments from two-dimensional Yang-Mills theory (YM2) and a duality 
symmetry found in the CSM, we find that the functional form of the measure must 
satisfy a host of properties. We find a function that satisfies these properties and 
conjecture that it is unique. 

Our starting point will be U (N) YM2 on the torus. It is by now well known 
that its partition function is given by the sum[8] 

Z = Y^e- LT ^' (2) 
R' 

where R' are the possible U(N) representations, L and T are the lengths of the 
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two circles comprising the torus, g is the YM2 coupling and C<ir* is the quadratic 
Casimir for representation R'. R' is described by a Young tableau with rij boxes 
in row i, which satisfy the constraints that ni > rij+i. Since we are considering 
U(N), the number of boxes in a row can be negative. The Casimir is then given 
by 

C, R , = 1 £ + JV-2i + l) = l£ («, + 2 - (3) 

Letting pj = rii + (iV — 2? + l)/2, then pi satisfies pi > p{ + \ + I. From this it is 
clear that the partition function in (2) is the same as that for iV nonrelativistic 
free fermions on a circle where the length of this circle, L', and the fermion mass 
m satisfy the relation mL' 2 = l/(g 2 L)[9}. T is the Euclidean time in this problem. 
Next insert a Wilson loop in representation R along the Euclidean time direction 
of the torus. This corresponds to putting a fixed color source at some point on the 
original YM2 circle. The loop modifies the partition function to 

Z R = J2 N R'R'R e ~ LT92C2R ' ( 4 ) 

R' 

where N^^^ is the number of singlets in the tensor product R\ x R2 x -R3. If R is 
the (I — 1)N symmetric representation of SU(N) then N R/ j^, R = 1 if pi > I +pi + \. 
Hence, the allowed states are precisely the states found in the CSM with coupling 
'[10]- 

The fact that the CSM is related to YM2 suggests that other techniques from 
YM2 might be applicable, in particular the string interpretation. This is equivalent 
to a formulation in terms of creation and annihilation operators. To achieve this, 
perform the standard redefinition of the wavefunction, \1/ = A'/, where A is the 
Vandermonde for Z{ = e 1 . The function / is a combination of YM2 characters, 
satisfying[2] 

H f = ^!T2UY. j^fM^ - z ^)f + 5> 2 ^ + *w ) = E f (5) 
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where d{ = d Zi . We can classify the states in terms of representations R' of U(N) 
(just as in YM2) as x(-R', I) with energy E(R',l). These can be thought as l- 
deformed characters and Casimirs. These are also referred to as Jack polynomials 
in the mathematical literature [11]. Note that 1) = XR is the usual U(N) 
character. / can be expanded in terms of the Schur functions W n = ^ zf, giving 

k 

f = ^2f{ni,...n k }Y[W ni = ^2f{m,...n k }\ni,...n k ) (6) 

k i=l k 

where f{rii} are expansion coefficients. Since here we will be interested only in 
chiral states, all rii can be chosen positive. This is the "string" picture, where 
the state is interpreted as a collection of strings with winding numbers rij. 

Defining string creation and annihilation operators \n) = a+|0) and n|0 >= a n \n) 
which obey the commutation relations [a n , = n8 m ^ n and interpreting the action 
of (5) on (6) in terms of a n and a+, we get [12, 13] 

H = -^jjy {(1 - l)na+a n + lNa+a n } + {ia+a+ a m+n + a+ +n a m a n } 

n>0 m,n>0 

(7) 

The terms inside the double sum in (7) are respectively a splitting term and a 
joining term. This string theory is nonlocal on the world sheet, due to the existence 
of the term na^a n , proportional to the square of the length of each winding string. 

A duality relation [11, 16] between the models with parameters I and l/l can 
be established using the Hamiltonian in (7). To this end, define the new "dual" 
operators a+ = — |a+ and a n = —la n satisfying the same commutation relations 
as the original ones. In terms of the new operators, the Hamiltonian becomes 

-I Tv^r,. 1, 1 



H = 



- j)na+a n + j{-lN)a+a n ] 



mL' 2 - 



I 

m,n>0 



So, apart from an overall constant —I, the Hamiltonian in terms of the new op- 
erators is identical to the old one but with / turned into l/l and iV turned into 



4 



—IN. Therefore, any solution of the l/l problem, expressed as a polynomial of 
creation operators acting on the vacuum, will also be a solution of the / problem 
upon substitution of the dual operators in its expression. But a+ creates the state 
—lW n . Further, turning W n into —W n in the expression for the character of a rep- 
resentation R' changes it into the character of the dual representation R', where 
the tableau of R' is that of R' with rows and columns interchanged. Hence, one 
has the duality relations 

X (R', Z; {W n }) = X (R', y , {-lW n }), E(R', I, N) = -IE(R' , j, -IN) (9) 

Changing the positive integer N into —IN may look unphysical; however, N ap- 
pears simply as a parameter in the Hamiltonian (its other manifestation is in 
nonperturbative effects). Therefore we can formally take it to be — IN in all ex- 
pressions explicitly involving it. For generic I, duality maps particles into holes in 
the pseudomomentum spectrum. Since n particles can annihilate with In holes, 
this justifies the mapping I — > l/l. Obviously, duality is broken by nonperturbative 
effects, since particle excitations can have arbitrarily high momentum while holes 
cannot be excited any further than the bottom of the Fermi sea. 

In order to calculate correlation functions we must know the inner product 
of the wavefunctions in terms of the correct measure. It can be shown that in 
the large-iV limit < n\m >= j5 m ,n- In fact, this relation implies that in this 
limit, a n = la^ (since < n\m >= < 0|(a+)ta+ 10 >). The Hamiltonian (7) is 
Hermitian under this measure, and therefore the eigenstates are orthogonal with 
respect to this measure. To obtain the exact measure with direct calculation is 
nontrivial for general I. The alternative is to determine it indirectly, through the 
jacobian of the change of variables from the old wavefunction ^{Oj) to the new one, 
viewed as a function of the Schur variables W n . The evaluation of this jacobian 
is achieved with standard collective field theory techniques Fourier transformed to 
momentum space[14,15,13]. We simply quote the result (appropriately modified 
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for the periodic Sutherland problem) 

OO 

, nJ = E; jj^, + a _ 1)E( _ 1)t i E Wm ^ Wn (10) 

n>0 ' fc=2 V ; V =0 

This result is exact perturbatively to all orders in 1/N, and it contains an infinite 
number of vertices, of increasing order in 1/N. The inner product of any two states 
expressed as functionals of W n , is given by integrating over all W n with the measure 
J j J J Y\n dW n . In particular, the overlap < n\ . . . n^n^i . . . nfc +m > is a k + m- 
point correlation function of an interacting one-dimensional field theory with action 
(10). Since (10) is invariant under W n — > —lW n , N — > — IN, I — > 1//, the norm 
and the inner product of the states is invariant under duality transformations. The 
leading quadratic part of (10), of order N°, reproduces the leading-order result for 
the overlap (n\m). Since, however, we are interested in the thermodynamic limit 
where the particle density remains finite we would like to compute the result to all 
orders in 1/N. 

YM2 provides a useful framework for determing selection rules in the density 
correlation functions. The density operator for a system of iV fermions is given by 

p(u) = (I/)" 1 *0 = (Lr 1 ££e- <m,/I/ e <nfl4 . (11) 

i n i 

But the sum over i is just the trace of an N by iV unitary matrix whose eigenvalues 
are e tdi . Hence p(u) is p(u) = (2nL')^ 1 J2 n e~ mu ^ L 'W n . The density correlation 
function is found by inserting onto the torus two spatial Wilson loops, p(u) and 
p{v), which are separated by Euclidean time T and which cross the Wilson loop, R. 
The loops divide the surface of the torus into two regions and to find the correlations 
on the ground-state the area of one of these regions is taken to infinity. Using 
arguments of Migdal and Rusakov[8], the density correlation function reduces to 

(p(u,t)p(v,0)) = —^Y,z~ m{U ~ V)IL ' j dUdVdWj^dR'dGXGiUWVW) 
XR'(UV^) x X r(W) tr(U n ) tr(V^)e- tLC « 

(12) 
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G is the so called "staircase" with I — 1 boxes per step and dji> is the dimension 
of representation R' . The integrals are over the U(N) group volume. The Schur 
function W n can be expressed as the sum of U (N) characters 

W n = X [n](U) ~ X[n-ld U ) + X[n-2,1,1] ( U ) + -(-^TX^.l] (U) (13) 

where the numbers in the square brackets correspond to the number of boxes in 
each row of the corresponding tableau. The only representations that appear in 
W n are the "corner" tableaux, that is, those that only have boxes in the first row 
or first column. 

In principle, the integration (12) can be carried out, but in practice it is highly 
tedious. The result of each integral will be given by some Clebsch-Gordan coeffi- 
cients. However, the contribution from R' is zero in (12) if the tensor product of G 
with all possible corner representations does not contain R' . This then leads to our 
selection rules, which are a generalization of those in [7]. Define the representation 
G p as the representation obtained by adding the boxes in row i of p to row i of G. 
Then it is straightforward to see that the allowed representations R' are given by 
G p , where p can have more than one box in each of the first / columns, but has at 
most one box after that. 

The correlation functions of greatest interest are those found in the large radius 
limit, where the overall fermion density is kept fixed. This implies that N/L' is 
fixed and therefore g 2 N 2 is fixed. This differs from the more traditional large 
N limit of Yang-Mills theory where g 2 N is fixed. With this new limit, a U(N) 
representation which has order N boxes in its tableau has a finite energy in the 
large iV-limit. It is useful to define new variables X{ = rii/N and y\ = fii/N, 
where rii is the number of boxes in the 2 th row and hi is the number in the i th 
column. To leading order in 1/N, the energy of a particular representation R' is 
E(xi,yi..yi) — 3 \ L ^1(^1 + — — Vi)j- From the selection rules, the 

density correlation function is given by I + 1 sums, where the first / are the sums 
over the number of boxes in the first I columns, and the other sum is over the 
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number of boxes in the first row. In the scaling limit, the sums become integrals 
and the integration limits are < x\ < oo and < y\ < y\-\ < y\-2--V\ < 1- It 
is easy to understand the nature of these sums by considering the Bethe ansatz 
picture for the CSM[2,17]. The ground state has I — 1 empty sites between each 
filled site. A representation R! corresponds to the situation where fermion 1 moves 
over n\ sites, fermion 2 moves over I sites, as do the next hi — 2 fermions. Then 
the next fermion moves over only I — 1 sites, and the logic is continued all the way 
down to fermion h\ which moves over 1 site. This is then equivalent to exciting 
one particle out of the fermi sea and leaving / holes, lyi — 1/2 is the momentum 
of i th hole, and x\ + 1/2 is the momentum of the particle. Hence the correlation 
function should be given by 

oo l yi-i 

(p(u,t)p(v,0)) = ( 2 ^y)2 Re J dXl J dyi " ' J ^^I'fi-SwO / 14 x 



iNix.+^V^iu-v)/^ ' -tN 2 [x 1 (x 1 +l)-l^.y i (l-y i )]/(2mL' 2 ) 



e 



It is now just a question of determining p(x\, yi-.yr, /)• To this end, we make the 
following observations: (i) The measure should be invariant under an interchange 
of the yi, since this just corresponds to interchanging the momenta of the holes. 
(ii) The measure should be invariant under x\ — > — I — xi, yi — > 1 — yi, since this 
corresponds to a parity transformation. (Hi) The integrand should vanish as yi — > 
yj, since there can not be two holes in the same place, (iv) Since (W- n W n ) = N 
if n > IN, the measure is nonzero and finite as x\ — > oo. 

For a particular eigenfunction of the Hamiltonian, the contribution to the mea- 
sure is proportional to {W n \^) {^l^) -1 (ijj\W n ) . There are some known eigenf unc- 
tions for the CSM, such as the n-antisymmetric wave-functions, which are given 
by 

n partitions k=l 

The sum is over all partitions of n. We can then make the additional observa- 
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tions: (v) Using the Jacobian in (10) the contribution of ip n to the correlation 
function is C(Z)n( 1-1 /0[i + (I - I)n/{IN)] + 0(l/iV 2 ), where C(Z) is an I dependent 
constant, (vi) By duality, the contribution of the 'symmetric' representation, if) n 
is C(l/Z)n( 1_/ )[1 + (/ - I)n/(IN)] + 0(l/iV 2 ). These six constraints are highly 
restrictive and the only function that we could find that satisfies them is 

u(x, „, v,-i) - c-' n^-te - w) 2/l fa(x 1 + ;))'-'fa + E,») 2 (15) 

/.(*i,»..»,o-c ^(l^),.!/. n ,:fe+w 2 (15) 

That this function satisfies the first four constraints is obvious. The last two 
require further explanation. The region of the integrand that comes from the 
antisymmetric wave function is 1/N ~ yi « y\ « 1 for i > 1, and 1/N ~ 
x\ << y\. Since we have replaced sums by integrals, the I + 1 integrals are 
integrated over a region of width 1/N. We then find that the contribution of ijj n 
is ~ (iV) 1 " 1 /^ 1 /' /(l —y) 1 ^ 1 / 1 , and the behavior agrees with (w). Likewise the 
contribution from vf) n comes from 1/N ~ ?/j << xi << 1, hence after integrating 
the variables over a width 1/N, one finds the behavior is ~ 
which agrees with (vi). 

In fact there are a few more checks that we can make. First, (14) and (15) 
agree with the previously known cases of I = 1 and / = 2 [4]. Second, note that the 
contribution from W n comes from the region r = x\ + J2iUii where n = tN. We 
can pick off this contribution by inserting 5(r — x\ — Hi) into the integrand. If 
r > /, we can do the x\ integration, leaving the integrals over yi saturated. But at 
equal time, we know that (W n W~ n ) = N for n > IN. Therefore, if (15) is correct, 
then it must be true that 



i/n (W ( 1 - W )A^.-^n(»*-^-^-^o l )'-M-c 

1 1< ~J 

(16) 

where Y = YliVi an< ^ C is a constant. The constant is determined by taking 
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r — > oo, in which case (16) reduces to the Selberg integral and C is given by 



We have been unable to prove (16) directly, but we computed the / = 3 case 
numerically and were able to show that it agreed with C up to six significant digits 
all the way down to r = 3.5. Finally, we observe that as x\ — » oo, the measure 
reduces to the measure recently found by Forrester for the two-particle Green 
function[18,6]. This is sensible, since sending the particle to infinite momentum 
should be equivalent to removing it altogether. 

Using the duality symmetry we can find the correlation function for the case 
1/7. In this case, we have I integrals over variables x% that run from to oo and one 
integral over y\ from to 1. The correlation function is then found by substituting 
—N/l for N, —lyi for x\ and — Ixi for m in the measure and energy. The result for 
/ = 1/2 agrees with the result found previously [4]. 

Note added After this work was completed, we learned that Ha has been able to 
prove these results using properties of the Jack polynomials [19]. After this paper 
was submitted we received another preprint proving these results [20] 
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